function [y] = inverse_au_carre(x1,x2)
    y = 1 ./ (x1+x2).^2;
end

